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ABSTRACT 

We investigate the radial transport of magnetic flux in a thin accretion disc, the turbu- 
lence being modelled by effective diffusion coefhcients (viscosity and resistivity). Both 
turbulent diffusion and advection by the accretion flow contribute to flux transport, 
and they are likely to act in opposition. We study the consequences of the vertical 
variation of the diffusion coefficients, due to a varying strength of the turbulence. For 
this purpose, we consider three different vertical profiles of these coefficients. The first 
one is aimed at mimicking the turbulent stress profile observed in numerical simula- 
tions of MHD turbulence in stratified discs. This enables us to confirm the robustness 
of the main result of Paper I obtained for uniform diffusion coefficients that, for weak 
magnetic fields, the contribution of the accretion flow to the transport velocity of mag- 
netic flux is much larger than the transport velocity of mass. We then consider the 
presence of a dead zone around the equatorial plane, where the physical resistivity is 
high while the turbulent viscosity is low. We find that it amplifies the previous effect: 
weak magnetic fields can be advected orders of magnitude faster than mass, for dead 
zones with a large vertical extension. The ratio of advection to diffusion, determining 
the maximum inclination of the field at the surface of the disc, is however not much 
affected. Finally, we study the effect of a non-turbulent layer at the surface of the 
disc, which has been suggested as a way to reduce the diffusion of the magnetic flux. 
We find that the reduction of the diffusion requires the conducting layer to extend 
below the height at which the magnetic pressure equals the thermal pressure. As a 
consequence, if the absence of turbulence is caused by the large-scale magnetic field, 
the highly conducting layer is inefficient at reducing the diffusion. 

Key words: accretion, accretion discs - magnetic fields - MHD - ISM: jets and 
outflows - galaxies: jets. 



1 INTRODUCTION 

Jets or outflows are observed in a variety of objects where 
an accretion disc rotates and slowly accretes onto a cen- 
tral object. This central object varies considerably in mass 
from a supermassive black hole in AGNs, to a stellar-mass 
black hole or a neutron star in X-ray binaries, a white 
dwarf in cataclysmic variables or a pre-main sequence star 
in T Tauri star systems. The presence of outflows in all 
these objects at very different scales suggests a deep con- 
nection between the accretion process through an accretion 
disc and the launching process of the outflow, which is fur- 
ther suggested by a correlation between the observational 
signatures of ejection and accretion. Theoretical models gen- 
erally appeal to a large-scale magnetic fleld to accelerate and 
coUimate the outflow into a jet. The magnetic fleld serves 



as an agent to extract the energy contained in the rota- 
tion o f the accretion disc in the magnetocentrifugal mecha- 
nism (iBlandford fc Pavnd|l982l ') or of the central black hole 



IBlandford fc ZnaieklllQTTi r It also has the advantage of pro- 



viding a natural collimation of the outflow. But can an ac- 
cretion disc possess a strong magnetic field as required by 
these models? And what is its origin? 



One possibility is a dy namo arising from the ma gnetoro- 
tational instability (MRI. iBalbus fc HawlevI l|l99ll ')'). How- 
ever, the dynamo process would need to be very efficient 
to produce the right strength of the magnetic field. It also 
remains unclear whether a large-scale magnetic field can be 
created in this way, because the correlation length of the tur- 
bulence is probably of the order of the vertical scale height 
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of the disc, which is muc h smaller than its radius (see the 
discussion in ISpruitll2010l ). 

In another scenario, the gas supplied to the disc contains 
an initially weak magnetic field which is transported inwards 
by the accretion flow. The magnetic fleld is strongly ampli- 
fied as magnetic flux accumulates in the central part of the 
accretion disc, and could thus reach the required strength 
to launch an outflow. The efficiency of this scenario depends 
upon the relative efficiency of the inward advection of mag- 
netic fiux by the accretion flow and of the diffusion of the 
magnetic fleld by the action of the efl'ective turbulent resis- 
tivity, which tends to counteract the advection process and 
transport magnetic flux outwards. Unfortunately, early the- 
oretical studies suggested that in thin accretion discs the 
diffusion process is much more efficient than the a dvect ion 
process ow et al. 1 11994 iHevvaerts et~al]|l996l ). Indeed 
the diffusion arises from the vertical diffusion of the radial 
component of the magnetic field across the accretion disc. 
In a thin disc, this takes place on a much smaller scale than 
the radius, which is the scale on which the effective viscosity 
transports angular momentum thus driving the advection. 
This scale difference makes the diffusion much more efficient 
than the diffusion. This raises doubts over the existence of 
a significant large-scale magnetic fleld in accretion discs and 
is therefore a major problem for models of jets. 

The aforementioned results were, however, obtained 
from a crude vertical averaging, and it was re- 
alised later that the vertical str ucture of the accre- 
tion dis c might ch ange the picture (Ogilvie & Livio 200l|; 
[Bisnov atvi-Kogan fc Lovelace 2007; GuUet & Ogilvie 201^. 
The purpose of this series of articles is to investigate the in- 
fluence of the vertical structure of the disc on the transport 
of magnetic flux, the effect of turbulence being modelled 
by effective diffu s ion co efflcients (viscosity and resistivity). 
iGuilet fc Qgilvid (|2012| ) (hereafter caUed Paper I) was re- 
stricted to the idealised case of uniform diffusion coefficients. 
It was found that even in that simple case the vertical pro- 
files of the velocity and magnetic field have a significant 
impact on the transport of magnetic fiux (particularly for 
magnetic fields that are weak in the sense that the mag- 
netic pressure is small compared to the thermal pressure 
in the equatorial plane). This article is devoted to study- 
ing the impact of a vertical variation of the diffusion coef- 
ficients. It is indeed expected that the turbulence strength 
varies with the distance from the midplane, as shown by nu- 
merical simulations of stratifi ed discs (e.g. iMiller fc Stonel 
I2OOOI : iFromang fc Nelsonll2006l '). This should result in a ver- 
tical variation of the effective diffusion coefficients. Three 
different physically motivated profiles are considered in this 
paper: 

• The first one, inspired by MHD simulations of fully tur- 
bulent discs, is mostly aimed at testing the robustness of the 
results obtained in Paper I with a more realistic description 
of the vertical structure of the turbulence. 

• Secondly, we consider the presence of a dead zone close 
to the equatorial plane, where the (physical) resistivity is 
large owing to low ionisation. As a consequence the MRI- 
induced turbulence is switched off, resulting in a small effec- 
tive viscosity. Dead zon es are expected to exist in some part s 
of protoplanetary discs (|Gammielll996l : lFromang et al.ll2002l ) 



and could therefore influence the transport of magnetic fleld 
in these objects. 

• Finally, a third profile of the diffusion coefficients 
is considered in order to te s t the scenario proposed by 
iBisnovatvi-Kogan fc Lovelace! (|2007l ) to reduce the diffusion 
of the magnetic flux. It relies on a non-turbulent layer at 
the surface of the disc, where the magnetic field is strong 
enough to stabilise the MRI. According to this scenario, 
the reduced turbulent diffusion in this layer would allow the 
magnetic field lines to bend without much diffusion. How- 
ever, our self-consistent calculation of the vertical structure 
with a non-turbulent surface layer shows that this layer is in- 
efficient at reducing the diffusion unless it extends to where 
the magnetic field is not expected to stabilise the MRI. 

In Section [2] we summarise the formalism developed in 
Paper I and used in the rest of the paper to compute the 
vertical structure of an accretion disc and the transport rate 
of magnetic flux. In Section |3l we consider a vertical profile 
of the turbulent viscosity that mimics the results of MHD 
simulations of fully turbulent discs. In Section |4l we model 
a non-turbulent dead zone around the equatorial plane. Fi- 
nally in Section \5\ we study the effect of a non-turbulent 
surface layer on the diffusion of the magnetic flux. We dis- 
cuss the results and conclude in Section]^ 



2 FORMALISM 

2.1 Assumptions and governing equations 

In Paper I, we derived a formalism enabling us to compute 
the transport of magnetic flux while taking into account the 
vertical structure of an accretion disc. In this section we 
summarise the approach and give the governing equations; 
the reader is referred to Paper I for more details on their 
derivation. The formalism is based on a systematic asymp- 
totic expansion in the limit of a thin accretion disc. For con- 
sistency in the expansion, as well as to obtain a linear system 
of equations which can be more easily solved, we further as- 
sume that the magnetic field is dominated by its vertical 
component. This amounts to neglecting the compression of 
the disc due to magnetic pressure. While this is an important 
limitation for strong magnetic fields, we showed in Paper I 
by comparing with direct numerical simulations that this ap- 
proximation gives good results for the range of magnetic field 
strengths considered in this paper. We assume a point-mass 
potential and therefore circular Keplerian orbital motion at 
leading order. An isothermal equation of state is assumed 
for simplicity. We look for solutions that are stationary on 
a dynamical timescale but may evolve on the longer viscous 
or resistive timescales. 

In order to model the effects of turbulence, we use ef- 
fective turbulent diffusion coefficients: a viscosity v and a 
magnetic diffusivity (proportional to resistivity) rj. The stan- 
dard alpha prescription is used for the viscosity: v = aCsH, 
where Cs is the sound speed and H is the vertical scale- 
height H = Cs/Ok with flK the Keplerian angular velocity. 
We allow for a vertical (but not radial) dependence of the 
a parameter through a — qo5(C)- The resistivity is then 
related to the viscosity through the magnetic Prandtl num- 
ber V = v/r], the vertical profile of which can also be freely 
chosen. The normalisation of the diffusion coefficients (i.e. 
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the value of qo) is determined by requiring that the MRI is 
made marginally stable by the diffusion. This ensures that 
the stationary solutions found are not unstable, and prevents 
unphysical effects such as multiple bending of the magnetic 
field lines that arise when unstable MRI modes exist (see 
Paper I for a discussion of this assumption). As already men- 
tioned, the vertical profile of the diffusion coefficients can be 
freely chosen, and three different profiles are considered in 
Sections El Hand El 

We use dimensionless variables defined as 



ters: 



pH 
E ' 

r Vr 

17 & ' 

H Cs ' 
r Br 

H ' 



(1) 
(2) 

(3) 

(4) 

(5) 



where p is the gas density, E the surface density, vk = i'Q.k 
the Keplerian azimuthal velocity, Vr and the radial and 
azimuthal components of the velocity, and Br and the 
radial and azimuthal components of the magnetic field. The 
vertical spatial coordinate z is normalised by the disc scale 
height, defining 



z 

H' 



(6) 



The vertical component of the magnetic field B^ does not 
depend on z in the limit of a thin disc and is therefore a free 
parameter of the problem. In the absence of an outflow due 
to the low inclination of the magnetic field lines, the vertical 
component of the velocity is negligibly small. 

With the isothermal equation of state, the density pro- 
file is a Gaussian function: 



P = 



I 



-CV2 



(7) 



Under these assumptions, the vertical structure of Ur, 
u^, hr and is determined by a quasi- local effectively sta- 
tionary problem, where global gradients appear as source 
terms. The governing differential system is 

I 13 

— -9c (pctdQUr) — 2u^ — -—d^br = - 



+ 



Db_ 

Pop' 



(8) 



1 113 

-de {padQU^) + -Ur - -r-zd^b^ = -a 
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+Dh 1 + 



dlnQ 
dlnC 



D^s - D„C 



-de {^9cbr 



d^Ur = —DbOq 



(9) 



(10) 



{^dcb^) - dcu^ + ^br = 0, (11) 
where we have defined the following dimensionless parame- 



/3o = 



Dh = 



Db = 



PO ECa 

Bi H 
dlnH 
dlnr ' 

2Dh - 

dlnB, 
dlnr 



3 ainS 
2 ainr ' 



(12) 
(13) 
(14) 
(15) 



Here /3o corresponds roughly to the midplane value of the 
plasma /3 parameter (the ratio of the thermal pressure to 
the magnetic pressure); more precisely, the two are re- 
lated by /3((" = 0) = \/2/tt /3o. The parameter D^e equals 
i91n(!/E)/91nr, given that u = aH^Q and a does not de- 
pend on r. 

The boundary conditions at — >■ ±00 are motivated by 
the absence of outflow, owing to the low inclination of the 
field lines. As a consequence the magnetic field at infinity is 
force-free, giving the following conditions at — >■ ±00: 



pUr — >■ 0, 



pu^ — > 0, 



br-{DBC±brs)^0, 



bet, - (i^'ps) 0. 



(16) 



(17) 



(18) 



(19) 



The inclusion of a non-vanishing azimuthal component of 
the magnetic field as — > ±00 allows for the possibility of a 
magnetic torque between the disc and the external medium. 
If the disc is isolated and has no outflow, then 6^a should 
be zero; however, the parameter fe^s allows us, if we wish, 
to mimic in a simple way the effect of angular momentum 
removal by an outflow or by the connection to an external 
medium. These boundary conditions are homogeneous, ex- 
cept for the linear source terms proportional to _D_b, bm, and 
b^B in equations (|18p and (|19p . 



2.2 Form of the solution 

The system of equations given above is a linear system of 
ordinary differential equations for itr(C), ^'<#>(C)i ^'■(C) a^nd 
fe0(C); radial derivatives appear only as parameters of this 
quasi-local problem. As a consequence of the linearity of 
the equations the solution depends linearly on the source 
terms, appearing either on the right-hand side of the sys- 
tem of equations or as a non-vanishing boundary condi- 
tion at infinity (6rs and fc^s)- Thus, the general solution 
X = {ur,u^,br,b^} is a linear combination of the solution 
vectors corresponding to each source term: 



X = Xk + XdhDh + X di^sDvj: + XdbDb 

-\-Xhrsbrs ~t" Xh(j)sb(f)s^ 



(20) 



where Xdh is the solution vector corresponding to the 
source term proportional to Dh and so on. Xk is the so- 
lution vector corresponding to the source terms that are 
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independent of Dh, -D^s, Db, brs, and b^s; these terms con- 
tain the radial derivative of the leading-order angular veloc- 
ity which is assumed to be Keplerian (hence the notation), 
as well as a term describing the vertical dependence of the 
gravitational potential and other geometrical terms coming 
from differentiating factors depending on r in the viscous 
term and the radial pressure gradient. Furthermore, we de- 
fine 



where 



' hyd 



(21) 



being the solution corresponding to the hydrodynamic 
source terms with the standard parameters Dh ~ 1 (rele- 
vant to a disc with a constant aspect ratio H/r) and Di/e ~ 
(relevant to a steady accretion disc far from the inner bound- 
ary). For these parameters, we thus have 



X 



hyd 



XdbDb + Xtrsbrs + X f, 



(22) 



The linearity of the problem makes the exploration of 
the parameter space and the physical understanding of the 
solutions much easier. Indeed, given the marginal stability 
hypothesis, the solution depends in a nonlinear way only on 
the two parameters /3o and V (which is assumed to be unity 
in this article, except inside the dead zone as described in 
Section |4]). For each pair of values of these two parameters, 
one needs to compute the six solution vectors Xk, Xdh, 
XdvT., Xdb, Xbrs and -X^b^s, each of which can be rep- 
resented by plotting the profiles of Ur, U4,, for and fo^. The 
general solution is then just a linear combination of these 
solution vectors with the appropriate coefficients. To solve 
the problem numerically, we use a spectral method based 
on a decomposition on a basis of Whittaker functions as 
described in Paper I. 



2.3 Transport velocities from vertical averaging 

Once the vertical profiles of the velocity and magnetic field 
have been computed, the transport rates of mass and mag- 
netic flux can be obtained and expressed in terms of a trans- 
port velocity. These transport rates determine the evolution 
of the disc on a viscous or resistive timescale. The dimen- 
sionless mass transport velocity is, from equation ([9]), 



3 4 
pUr AC, = --a(l + 2D^T.) + -rb-tys, 
^ PO 



(23) 



where a is the density-weighted vertically averaged dimen- 
sionless viscosity: 



pa(C)dC. 



(24) 



For a vanishing fo^s this is consistent with the familiar ex- 
pression from the standard theory of a Keplerian accretion 
disc, 



(25) 



where the overbar refers to a density-weighted vertical av- 
erage. 

The integrated version of the radial component H10|) of 
the induction equation is 



= Ur + Tpidc^br — D b) = COUSt, 



(26) 



r 



(27) 



is the dimensionless magnetic flux transport velocity. It can 
be evaluated using equation (I26p at any convenient height. 
In particular, the force-free conditions at ^ — >■ ±00 ensure 
that the radial velocity at this location equals u^. This is 
useful to remember for the visual inspection of the verti- 
cal profiles displayed in Figures [3l ID |9] [10] and [121 A useful 
way of vertically averaging the induction equation to obtain 
the magnetic flux transport velocity is the following. Divid- 
ing equation (|26[) by the dimensionless resistivity a/V, and 
integrating between ^ — and we obtain: 



-Di 



+ Ur 



(28) 



where fj is an average dimensionless resistivity defined in 
the following way (it is actually the inverse of the height- 
averaged conductivity) : 



(29) 



and Ur is the average of the radial velocity weighted by the 
conductivity (l/j?): 



1 



C -p 
— Ur aQ . 



(30) 



This averaging procedure shows that the contribution of 
advection to the velocity at which the magnetic flux is trans- 
ported is the vertically averaged radial velocity weighted by 
the electrical conductivity. This average may be very dif- 
ferent from the density-weighted aver age Um that describes 
mass transport (jOgilvie fc Lividl200ll ). 

Note that equation (|28|l is valid for any height C, up 
to which the averaging procedure is performed, and which 
can be chosen in the most convenient way. Integrating up 
to an infinite height may not be very useful as the force-free 
regime would dominate the average and one would simply 
recover that the radial velocity at infinity equals u^. We 
showed in Paper I that it is more meaningful to perform the 
average up to the height C,b, defined as the height at which 
the magnetic and thermal pressures are equal: 



CB = ^/ln(-A? 



(31) 



Indeed, we showed in Paper I that for large enough val- 
ues of /3o (> 10'^) the vertical profiles can be described by a 
two-zone model, where the magnetic field is approximately 
passive for C, < C,b and force- free for C, > C,b- The velocity 
&t C, < C,B is not influenced by the magnetic held, implying 
that Ur can be obtained from the purely hydrodynamical 
problem. The value of for at = C,b is a bit less straightfor- 
ward to obtain because br can have a jump at the transition 
with the force-free region, where the boundary conditions 
at infinity can be applied. In Paper I, we found that for the 
source terms governing the diffusion of the magnetic field 
(fors and Db) the value of for at ^ = C^g was a factor 3/2 
larger than that at dictated by the boundary conditions: 
br{(,B) ~ (3/2)(fors + DbCb)- We find in the following sec- 
tions that this numerical factor can depend on the vertical 
profile of the diffusion coefficients but remains of order unity. 
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Figure 1. Model of a fully turbulent disc: vertical profiles of the dimensionless viscosity parameter a (left panel), and of the viscous 
stress {r(j> component of the viscous stress tensor: (3/2)pa, right panel) for a magnetic field strength /3o = 10^. The red line shows the 
profile used in Section [3] (equation |32]l, while the black line shows the case of uniform diffusion coefficients studied in Paper I. The vertical 
dashed lines show the height <^c below which the viscous stress is approximately uniform, and the dotted lines the height (^g above which 
the magnetic field is dominant. 




Figure 2. Diffusion coefBcients determined by the marginal sta- 
bility hypothesis as a function of the magnetic field strength. The 
black line correspond to a uniform a as studied in Paper I. Red 
lines correspond to the profile described by equation 1321 1 : the full 
lines show the minimum and maximum values of a, the dashed 
line shows the density-weighted average defined by equation 1241 1 
(relevant to the advection of mass), and the dash-dotted line 
shows the average resistivity defined by equation 1291 1 with the 
average being performed up to C = Cs (relevant to the diffusion 
of the magnetic field). 



3 MODEL OF A FULLY TURBULENT DISC 

As a consequence of the vertical structure of the density, 
the properties of the turbulence are not expected to be uni- 
formly distributed in the vertical direction inside an accre- 
tion disc. A large number of numerical simulations of 3D 
MHD turbulent discs with vertical stratification have been 
performed with both local and global models and can serve 
as a guide to build a more realistic model than the uni- 
form diffusion coefficients model of Paper I. Most of these 
simulations use a magnetic field configuration with zero net 
vertical flux because a significant net fiux tends to cause nu- 
merical difficulties with the vertical boundaries in the pres- 



ence of vertical stratification. One might, however, expect 
the stress profile in the case of a weak vertical field to be 
similar to that of these zero-net-flux calculations (at least 
close to the midplane, as long as /3 2> 1). Numerical simula- 
tions of stratifled discs with zero net vertical magnetic flux 
show that the disc is divided into two regions: the weakly 
magnetised main body of the disc at \z\ < 2 — 3H, and a 
strongly magnetised corona at \z\ > 2 — 3H. In the main 
body of the disc the turbulent stress is roughly independent 
of the distance from the midplane (e.g. ^Mil ler fc Stonell2000l : 
iFromang et al.l i2006'). In the strongly magnetised corona th e 
stress then decreases like the density (|Miller fc StondboOOl ) . 
Interpreting this vertical profile of the stress as a vertical 
dependence of the a viscosity leads to a profile of a that 
is inversely proportional to the density for |("| < 2 — 3 and 
constant above |c^| — 2 — 3. Such a vertical profile of a is 
probably more realistic than the uniform one for /3o 3> 1. 
We choose the following functional dependence of a: 



q(C) = Qo 



l + (e«/2_i)e-CV2' 



(32) 



where ao is the value of a at the midplane, and the pa- 
rameter ("c determines the height of the magnetic corona 
where a transitions from a profile approximately propor- 
tional to e'' (for < ^c) to a constant value (a ~ aoe''"^'^ 
for > (^c). In the following we use the value = 2, 
corresponding to a vertical variation of a by a factor of 



7.4. 



Figure [T] shows the vertical profile of the a parameter 
(left panel) as well the viscous stress (right panel) obtained 
for our fiducial value of the magnetic field strength /?o = 10*. 
The a parameter ranges from 0.12 at the midplane to 0.9 in 
the corona. The density-weighted average of a (relevant to 
advection of mass) is a = 0.2, while the average resistivity 
as defined in equation (|29[) (relevant to the diffusion of the 
magnetic flux) is fj = 0.31. The dependence of these various 
measures of the diffusion with respect to the magnetic fleld 
strength is shown in Figure [21 Note that for > 10'' - 10^ 
the value of a in the corona is the same as that obtained in 
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the case of uniform diffusion coefficients (and ranges from 0.5 
for weak fields to 1). This is explained by the fact that in the 
marginal stability hypothesis the value of a is determined by 
the largest-scale MRI channel modes, localised around ^ — 
^B, which is in the corona since ("_g > fc- As a result of the 
profile considered, the midplane value as well as the vertical 
averages of the diffusion coefflcients are significantly smaller 
(and therefore more realistic) than the value in the corona: 
for the weakest magnetic fields considered the midplane a 
becomes as low as 7 x 10~^, though this is still a bit higher 
than the value generally obtained in numerical simulations. 

Figure [3] shows the vertical profiles of the velocity and 
magnetic field corresponding to the source term brs, which 
describes the vertical diffusion of the radial component of 
the magnetic field through the disc. The solution obtained 
with the diffusion-coefficient profiles given by equation (|32p 
(shown with red lines) is compared to the solution obtained 
with uniform diffusion coefficients (shown with black lines). 
With uniform diffusion coefficients, the radial component of 
the magnetic field has a linear dependence on for \(\ < C^b- 
The shape is more complicated with the new a profile: the 
bending of the field line occurs preferentially close to the 
midplane (|("| < (c) where a is low. This can be understood 
using the fact that the magnetic fiux transport velocity 
expressed by equation (|26[1 is independent of the vertical 
coordinate. Indeed, for |(^| < (b the magnetic field behaves 
approximately passively, therefore the radial velocity van- 
ishes. The derivative of the radial magnetic field can then 
be directly expressed as 

V 

dthr = —u^ + Db, (33) 
a 

which shows that the bending occurs preferentially where 
the resistivity a/V is low. Note that this property is the 
physical reason behind the definition of the average resis- 
tivity fj to be the inverse of the average conductivity: more 
weight is given to the region where the conductivity is high 
because the magnetic field line bends preferentially there. 

Another difference between the two diffusion coefficient 
profiles is the amplitude of the jump in radial magnetic field 
around i^s: the jump is more pronounced with the new pro- 
file with Abr — 1, while it is Abr — 0.5 for uniform diffusion 
coefficients (see Paper I). This is intimately related to the 
larger bump in radial velocity around (b due to the require- 
ment that magnetic flux transport velocity be independent 
of (equation I26p . A description of the transition region 
(between the passive magnetic field regime at |(^| < (^b and 
the force-free magnetic field regime at |("| > ("s) was given 
in Appendices A and B of Paper I. It shows that the jump 
in br around (^b is related to the jump in azimuthal velocity 
through Abr = — 1 4—^ . The numerical solutions show that 
the azimuthal velocity jump is approximately twice larger 
with the new viscosity profile, which is in agreement with 
this relationship. The jump in azimuthal velocity is in turn 
determined by the requirement that the radial electric cur- 
rent be continuous across ("s: Au^ = {a/P)d(;btf,{<^g) (equa- 
tion 100 of Paper I). The larger jumps of br and in the 
new profiles can then be traced back to the fact that the 
azimuthal magnetic field is larger because the smaller resis- 
tivity near the midplane allows more shearing of the radial 
field into an azimuthal field. 

Figure [4] shows the vertical profiles of the velocity and 



magnetic field for the hydrodynamical source terms, describ- 
ing the advection due to angular momentum transport by 
the effective viscosity. At |("| < (c, the radial velocity is more 
negative in the case of uniform diffusion coefficients, which 
can be interpreted as due to the larger viscosity there com- 
pared to the more realistic profile. Note that for the latter 
profile, the radial velocity is even positive at \(\ < 0.78 with 
a small maximum value of Ur = 0.2. Such a positive ra- 
dial velocity near the midplane of an accretion disc where 
the averaged radial velocity is negative is sometimes called 
meridional circulation. It has already been found in viscous 
models of accretio n di scs with uniform a fo r example by 
iKlev fc LinI l|l992l l and lTakeuchi fc LinI l|2002l ). Such a phe- 
nomenon is absent from the uniform a case because of the 
large value of a (equation 87 of Paper I shows that it would 
appear at small q). At \(\ > ^c, the radial velocities of the 
two viscosity profiles become comparable as the viscosity 
also becomes comparable (see Figure [T]). As with the source 
term brs, the smaller resistivity near the midplane leads to 
an easier bending of the field lines and therefore larger val- 
ues of the radial and azimuthal magnetic field components 
(and therefore also a larger jump of the azimuthal velocity 
around (b)- 

Figure [5] shows the transport velocity of the magnetic 
flux for the different source terms (red lines) and compares it 
to the case of a uniform a (black lines). The transport veloc- 
ity due to the inclination of the field lines (brs) is normalised 
by the averaged resistivity t] as suggested by the averaging 
procedure described in Section [2.31 (it is shown in the upper 
left panel). Normalised in this way, the diffusion velocity is 
a bit larger with the new profile of a because of the larger 
jump in br discussed above. All in all this diffusion velocity 
is, however, smaller with the new profile due to the lower 
value of fj, as can be seen on Figure [3] by comparing the 
radial velocity at large ^. The diffusion velocity associated 
with the source term Db (radial gradient of magnetic field 
strength), normalised in the same way, is also larger with 
the new a profile (upper right panel). This is due to the 
same phenomenon: the jump in br is larger because the re- 
sistivity is reduced in the midplan^H- The transport velocity 
of the magnetic fiux due the hydrodynamical source term, 
shown in the lower left panel, is normalised by the trans- 
port velocity of mass. Thes e two transport velo cities were 
assumed to be the same in iLubow et ai] l|l994 ) and other 
papers, but we showed in Paper I that they can differ signif- 
icantly by a factor of up to 10 for the magnetic field strength 
considered. This is explained by the fact that the transport 
velocity of mass is dominated by the region close to the 
midplane (since it is a density-weighted average) while the 
transport velocity of magnetic fiux is significantly affected 
by the low-density region away from the midplane (since it is 
a conductivity weighted average) where the velocity is large. 
The same qualitative feature is found with the more realistic 
viscosity profile, but the ratio u^/um is even higher, reach- 
ing up to 25 for Po = 10* . This is a consequence of the larger 



The difference in diffusion velocity between the two a profiles 
is more pronounced for the source term Db than for the source 
term brs. This is because all of the diffusion is due to the jump 
in br in the case of Db , while part or most of it is not due to this 
jump in the case of brs- 
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Figure 3. Diffusion of a radially inclined magnetic field with a vertical profile of the viscosity corresponding to a uniform stress in the 
main body of the disc. The four panels correspond to the vertical profiles of the radial (left) and azimuthal (right) magnetic field (top) 
and velocity (bottom) when a radial magnetic field bra = 1 (—1) is imposed at the upper (lower) boundary (i.e. the solution vector 
^brs)- The red line corresponds to the diffusion coefficient profile of Section l3l (equation I32|l . while the black line is the uniform diffusion 
coefficient case studied in Paper I. The strength of the vertical magnetic field is /3o = 10^. The vertical dashed and dotted lines have the 
same meaning as in Figure [T] 




Figure 4. Same as Figure [3l but for the 'hydrodynamic' source terms (Dh = 1 and D^e = 0). 



difference between the (smaller) midplane velocity and the 
(similar) velocity in the corona. Finally the transport veloc- 
ity of magnetic flux due to an outflow (source term 6,^s) is 
shown in the lower right panel. It is found to be smaller with 
the new diffusion coeflicient profile by a factor of ~ 3. 

An important consequence of the relative efficiency of 
the diffusion and advection processes is the radial inclina- 



tion with respect to the vertical direction that magnetic field 
lines can sustain in a stationary situation. When the advec- 
tion process is due to the angular momentum transport by 
the effective viscosity, the inclination is expressed as (see 
Section 4.6 of Paper I) 

Brs ^ H Mhyd ,g^s 
-Bz r Ubrs ' 
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Figure 5. Contributions to the transport velocity of the magnetic flux, u^, as a function of the magnetic field strength, for different 
source terms. Upper left panel: diffusion due to the bending across the disc {bra source term). Upper right panel; diffusion due to the 
radial gradient of {Dg). These two transport velocities are normalised by the average resistivity defined in equation II29I I. Lower left 
panel : advection due to the turbulent viscosity ('hydrodynamic' source terms, Dh = 1 and -D^e = 0) normalised by the advection 
velocity of mass. Lower right panel: advection due to angular momentum removal by an outflow (fe^a). In all panels, the dotted black 
line corresponds to the case of uniform diffusion coefficients studied in Paper I, the red line corresponds to the viscosity profile giving a 
uniform stress in the main body of the disc (Section [Sjl and the blue line corresponds to diffusion-coefficient profiles modelling a dead 
zone (Section IJll. 



!> 40 





Figure 6. Radial inclination of the field lines at the surface of the disc obtained in a stationary situation. The left panel shows the 
inclination (Brs/Bz) in units of H/r such that the diffusion compensates the advection due to the turbulent viscosity ('hydrodynamic' 
source terms). The right panel shows the ratio —bre/b^s such that diffusion compensates advection due to an outfiow. In both panels, the 
dotted black line corresponds to the case of uniform diffusion coefficients studied in Paper I, the full red line corresponds to the profile 
with a constant stress in the midplane (Section O , and the full blue line corresponds to the dead-zone model (Section [4t. 



This inclination in units of H/r is shown in the left panel 
of Figure [6l Interestingly this inclination has a very similar 
behaviour with the new diffusion-coefficient profile as with 
uniform diffusion coefficients (the two differ by about 20%): 
it increases as the magnetic field strength is decreased, and 
reaches quite large values o f up to 50 (to be compared to 3/2 
used bv lLubow et al.lll994 ). This shows that this important 



result of Paper I is robust and remains true for more realistic 
treatment of the vertical structure of the turbulence. 

If the advection of magnetic flux is due to to angu- 
lar momentum loss in an outflow, the ratio of radial to az- 
imuthal magnetic field at the surface of the disc in a sta- 
tionary situation is determined by 

Bra Ubc^s /qc:\ 

R~ ^ '^7, — ■ ^ ' 
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Figure 8. Diffusion coefficients determined by the marginal sta- 
bility hypothesis as a function of the magnetic field strength for 
the dead zone model of Section |4] The black line shows the uni- 
form a case for comparison. Blue lines correspond to the profile 
modelling a dead zone given in equations I I36I I- II39I I: the mini- 
mum and maximum values of a are represented by full lines, the 
density-weighted average by the dashed line, and the average re- 
sistivity (Equation I29I I by the dash-dotted line. 



This quantity is shown in the right panel of Figure [S] It 
is smaUer with the new profile because of the lesser effi- 
ciency with which the accretion flow driven by an outflow 
can transport magnetic flux. 



4 DISC CONTAINING A DEAD ZONE 

The diffusion-coefficient profile studied in Section [3] models 
an accretion disc where MHD turbulence operates through- 
out the disc. This is expected for discs which are sufficiently 
ionised to enable the MRI to operate, such as the hot discs 
around black holes or a neutron star in AGN and X-ray 
binaries. The cooler discs around a T Tauri star may not 
be sufficiently ionised throughout their volume to sustain 
MHD turbulenc e, and a lamina r region called a 'dead zone' 
may then exist l|Gammielll996l ). Determining the transport 
of magnetic flux in such discs then requires a speciflc mod- 
elling of this dead zone, which is the purpose of this section. 
Understanding the strength of the magnetic field in proto- 
planetary discs is important f or the launching of o utfiows as 
outlined in the introduction (jPerreira et ahlliood ). but also 
for the theory of planet formation. Indeed dead zones are a 
probable formation site for planets, and recent studies show 
that the magnetic field could strongly affect th e migration of 
planets (|Terauemll2003l : iFromang et al. I l2005l . Guilet et al., 
submitted to MNRAS). 

A disc can sustain MHD turbulence if it has a suffi- 
ciently low physical resistivity, which is governed mostly by 
the ionisation fraction of the gas. Different physical processes 
are responsible for the ionisation depending on the region 
of the disc considered. In the inner parts of protoplanetary 
discs (r < 1 AU) the temperature can be high enough for 
the thermal ionisation of metals to provide sufficient ionisa- 
tion (if T > 1000 K). At larger radii, non-thermal processes 
are the main agent of ionisation: X-rays or ultraviolet fight 
emitted by the central star or possibly cosmic rays. These 



ionising radiations are able to penetrate only a finite surface 
density into the disc (of the order lOOgcm"^, though this 
depends on the energy of the photons or cosmic rays). As 
a consequence, they can ionise the whole vertical extent of 
the disc only at large disc radii, typically r > 5 — 10 AU. At 
intermediate radii, only the surface layer of the disc is suf- 
ficiently ionised to sustain MHD turbulence, while a region 
around the midplane has a too high resistivity and therefore 
remains laminar (the dead zone). Note that the ionisation 
fraction is very sensitive to the amount of small dust parti- 
cles and metal in the gas eous phase, and as a consequence 
is very model-dependent l|Fromang et al.|[2002l ). The extent 
of the dead zone is therefore also very uncertain and the 
numbers quoted above should be considered as only rough 
estimates. 

In order to model the vertical structure of the region of 
a protoplanetary disc that contains a dead zone, we choose 
the following functional form of the diffusion coefficients: 



a(lCl > Cent) 
a/niCl > Cent) 

a(ICI < Cent) 
a/ViK] < Cent) 



= ao 



ao/Vo, 

^t^cnt + (1 — i^cut)e~ 
ao/Vo 



(36) 
(37) 



[r?cut + (l-r,cnt)e-(ICI-C-)ViLt] 



(39) 



where ^cut is the height up to which the dead zone extends, 
i^cut is the factor by which the viscosity is reduced inside 
the dead zone, l/rjcut the factor by which the resistivity is 
increased, and Lcut is the width of the transition between 
the active layer and the dead zone. As an illustrative exam- 
ple, we choose the following values: (^cut = 2, fcut — W~^, 
Vcut ~ 10~^ and icut = 0.5. The viscosity and resistiv- 
ity profiles are shown in the left panel of Figure [7] for a 
magnetic field strength of /3o = 10'*. In the active layers 
above i^cut , the viscosity and resistivity are very close to the 
case of uniform diffusion coefficients. Inside the dead zone 
the resistivity increases by a factor 10 while the viscosity 
is decreased by a factor of 100. The corresponding viscous 
stress profile is shown in the right panel of Figure [T] In- 
side the dead zone it decreases to a minimum value which 
is about 10% of the maximum stress reached in the active 
layer. This is consistent with the result of numerical simu- 
lations of local disc models containing a dead zone. These 
find that the stress in the dead zone is dominated by the 
reynolds stress due to waves excited at the boundary with 
the act ive layers, and is about 10% of the stress in the active 
layers dFleming fc Stonell2003l ). 

The dependence with respect to the magnetic field 
strength of various measures of the diffusion is shown in 
Figure [8l The density- weighted averaged viscosity (relevant 
to the transport of mass) is reduced by a factor ~ 7 due to 
the low viscosity in the dead zone. The averaged resistivity 
is only slightly larger than the resistivity in the active layers, 
although it is much larger in the dead zone. This is a conse- 
quence of the averaging procedure used: the volume average 
of the resistivity would be much larger than this inverse of 
the averaged conductivity. 

The vertical profiles of the velocity and the magnetic 
field obtained for the source term brs, describing the verti- 
cal diffusion of a radial magnetic field, are shown in Figure|9l 
The dead zone model (blue lines) is compared to the uni- 
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Figure 7. Disc containing a dead zone: vertical profiles of the diffusion coefficients (left panel), and of the viscous stress {rij> component 
of the viscous stress tensor) (3/2)pa (right panel) for a magnetic field strength /3o = 10*. In the left panel, the full line shows the 
dimensionless viscosity a, and the dashed line the dimensionless resistivity a/V. In both panels, the blue lines show the profiles used in 
Section[4]to model a dead zone (equations | |36I1 - ||39| |). while the black line shows uniform diffusion coefficient case studied in Paper I. 
The vertical dashed lines show the height i^cut delimiting the dead zone, and the dotted line the height (^g above which the magnetic 
field is dominant. 
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Figure 9. Same as Figure [3] but with diffusion coefficients profiles modelling a dead zone (shown with blue lines). The vertical dashed 
lines show the vertical extent of the dead zone. 



form difTusion coefficient profile of Paper I (black lines). The 
main difference is that the large resistivity inside the dead 
zone prevents the magnetic field lines from bending there, 
such that the radial and azimuthal components of the mag- 
netic field remain close to zero for C, < Ccut. All the bending 
therefore occurs in the active layers at (cut < C < Cs) where 
the profiles are qualitatively similar to the uniform diffu- 
sion coefficient case; in particular the radial component of 
the magnetic field has a linear dependence on Cj" because the 
diffusion coefficients are uniform in the active layers. The 
transport velocity of the magnetic fiux (visualised as the ra- 
dial velocity at large |^|) is slightly larger in the presence 
of the dead zone, consistent with the slightly larger average 



resistivity. This diffusion velocity normalised by the aver- 
age resistivity is shown as a function of the magnetic field 
strength in the upper left panel of Figure [S] It is slightly 
smaller in the dead zone model than in the uniform diffu- 
sion coefficient case, due to the fact that the bump in br 
around ("b is slightly less pronounced in the presence of a 
dead zone. 

The vertical profiles of the velocity and magnetic field 
for the hydrodynamical source terms are shown in Figure fTOl 
Again the radial and azimuthal magnetic field components 
remain close to zero inside the dead zone due to the large re- 
sistivity there. The radial velocity in the dead zone is slightly 
positive for \C^\ < 1.2, and becomes comparable to the uni- 
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Figure 10. Same as Figure|4]but with diffusion coefficients profiles modelling a dead zone (shown with blue lines). The vertical dashed 
lines show the vertical extent of the dead zone. 



form diffusion coefficient case (tiiough slightly less negative) 
outside of the dead zone. The advection velocity of the mag- 
netic flux of the two models is similar. This can be inter- 
preted by the fact that the radial velocity in the dead zone 
(where the two models differ most) barely affects the con- 
ductivity weighted averaged radial velocity because the con- 
ductivity is very small. The ratio of the magnetic flux and 
mass transport velocities is shown in the lower left panel of 
Figured it is a factor ~ 7 larger in the presence of a dead 
zone, because the mass transport velocity is 7 times smaller 
than in the case of uniform diffusion coefficients while the 
magnetic flux transport velocity is similar in both models. 

The left panel of Figure [6] shows the radial inclination 
of the magnetic field lines at the surface of the disc in a 
stationary situation where the diffusion compensates exactly 
the advection. It is very similar in the dead zone model and 
in the uniform diffusion coefficient model. Indeed we have 
seen that the advection velocity is very similar, and that the 
diffusion velocity is only very slightly larger. 

The transport velocity of magnetic ffux due to an out- 
ffow is shown in the lower left panel of Figure [S] It is a factor 
~ 1.5 — 2 larger in the presence of a dead zone than with 
uniform diffusion coefficients. As a consequence the ratio of 
radial to azimuthal magnetic field when this transport com- 
pensates diffusion is also a factor ~ 1.5 — 2 larger (left panel 
of Figure [6| . 



5 NON-TURBULENT SURFACE LAYER 



In a s eries of four pa pers 'Bisnovatvi-Kogan & Lovelace" 
ll2007h.lRothstein fc Lovelace (2008 '). Lovelace et al. (2009), 
and Bisnovatvi-Kogan fc Lovelace! (j 20 121 ') have proposed a 
model in which the diffusion of the magnetic field is re- 
duced owing to the existence of a non-turbulent layer near 
the surface of the accretion disc. In this model all the bend- 



ing of the field lines takes place in this surface layer with- 
out causing any diffusion because the conductivity is very 
high (due to the lack of turbulence). This model however 
still lacks a proper calculation of the vertical profile of the 
magnetic field showing such an effect. Indeed the first two 
papers use a heuristic approach and do not attempt such 
a calculation. The last two paper s (|Lovelace et al.l l2009l : 
iBisnovatvi-Kogan fc Lovelacel |2012| ) describe a formalism 
that in principle takes into account a non-turbulent surface 
layer. But solutions are then only calculated in the particu- 
lar case where the diffusivity and density are uniform inside 
the disc, and where the diffusivity vanishes only in the wind 
outside the disc, which is not computed explicitly. As a re- 
sult there is no change of the radial magnetic field in the 
surface layer, and therefore the diffusive transport of the 
magnetic ffux is not decreased. In these papers the main in- 
gredient allowing a significant bending of the magnetic field 
lines inside the disc is an increased advection due to the 
presence of an outffow, rather than a decreased diffusion 
due to a non-turbulent surface layer. This possibility was 
discussed in paper I. The purpose of this section is t o test 
the idea proposed in Bisnovatvi-Kogan fc Lovelace! (|2007l ) 



and iRothstein fc Lovelac? ! 20081 ) that a non-turbulent sur- 
face layer can reduce the diffusion of the magnetic ffux. To 
this end, we perform a calculation of the vertical profiles of 
velocity and magnetic field in a disc with a non-turbulent 
surface layer and determine the associated diffusion. 

In order to model the non-turbulent surface layer, we 
choose the following vertical profile of a: 



a{\C\ < Ccut) = "0, 

adCI > Ccut) = "0 I^Qcut + (1 — Q!cut)e~ 



(40) 

(ICl-Ccut)Vi^l^) 



where Ccut is the height above which the viscosity and re- 
sistivity are reduced (and below which they are uniform). 
Above Ccut, a transitions from oq to Qcutao over a height 
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Figure 11. Diffusion coefficients in the model of a non-turbulent surface layer. Left panel: vertical profiles of the dimensionless viscosity 
parameter a for /3o = 10^. The red line corresponds to a position of the non-turbulent layer of (^cut = 4, the blue line to fcut = 3, while 
the black line shows the uniform diffusion case of Paper I. The vertical dotted lines show the height (^g above which the magnetic field 
is dominant. Right panel: dimensionless viscosity in the midplane ao as a function of the height fcut of the non-turbulent surface layer 
normalised by the height (^g at which the magnetic field becomes dominant. 



of Lcut- The magnetic Prandtl number V is assumed to be 
uniform and equal to 1. In the following we used the val- 
ues Lcut = 0.5, and Qcut = 0.1 (we checked that decreasing 
this factor further does not change qualitatively the results). 
The location of the transition to a non-turbulent state ("cut 
was varied to determine where it should lie for the non- 
turbulent layer to effectiv ely reduce the diffusio n of m ag- 
netic flux. As discussed bv lRothstein fc Lovelac3 (|2008l ) the 
physically motivated position of the non-turbulent layer is 
Ccut — Cb, where the magnetic pressure equals the thermal 
pressure. Indeed, above (s the magnetic pressure exceeds 
the thermal pressure and as a consequence the MRI is sta- 
bilised, resulting in a laminar state, although the reality may 
be more complicated. 

The left panel of Figure [TT] shows the vertical profiles 
of the diffusion coefficients for /3o = 10* and for two posi- 
tions of the non-turbulent surface layer: ^cut = 3 (blue line) 
and (cut = 4 (red line) which is close to the expected value 
of (b — 4.3. These two profiles are compared to the uni- 
form diffusion coefficient case shown with a black line. For 
Ccut = 4, the marginal stability analysis gives a very simi- 
lar diffusion coefficients value inside the disc for both of the 
uniform coefficient case and the non-turbulent layer model. 
This is consistent with the idea that above (b the MRI is sta- 
bilised because of the strength of the magnetic field: no dif- 
fusion is therefore needed in this region to achieve marginal 
stability. For (cut = 3 (< (s), the marginal stability analysis 
gives much larger diffusion coefficients inside the disc. This 
shows that for this small value of ("cut < Cs, the magnetic 
field in the region of decreased diffusion is not large enough 
to stabilise the MRI, such that the diffusion has to be in- 
creased to achieve marginal stability. This idea is further 
demonstrated in the right panel of Figure 1111 showing ao 
(the value of a inside the disc) as a function of the position 
of the non-turbulent layer (cut normalised by its physically 
expected value Cs- For Ccut > Csi c^o is unchanged by the 
inclusion of the non-turbulent layer. On the contrary, for 
Ccut < Cs, Qo increases strongly when Ccut decreases, show- 
ing that we are trying to decrease the diffusion in a region 
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Figure 13. Diffusion of an inclined magnetic field in the presence 
of a non-turbulent surface layer: transport velocity normalised by 
a in the main body of the disc as a function of the height at which 
the non turbulent layer starts normalised by (^b (the height above 
which the magnetic field is dominant). The non-turbulent surface 
layer decreases the diffusion only if it extends below (b , which is 
not expected to be true. 



where the MRI is not stable. As a conclusion, the marginal 
stability analysis confirms that the expected position of the 
non-turbulent layer is Ccut — Cs- 

Figure [T^] shows the vertical profiles of the velocity 
and magnetic field for the source term brs, correspond- 
ing to the diffusion due to the radial inclination of the 
field lines at the surface of the accretion disc. The so- 
lution with the physically expected position of the non- 
turbulent layer is shown with red lines. It lies very close 
to the profiles corresponding to the uniform diffusion coeffi- 
cient case (black lines). This shows that the non-turbulent 
surface layer at its expected position barely affects the dif- 
fusion of the magnetic field, thus con tradicting the model of 
iBisnovatvi-Kogan fc Lovelace! ||2007|). Results simila r to the 
expectations of Bisnovatvi-Kogan fc Lovelace! (|2007! ) are ob- 
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Figure 12. Magnetic field diffusion in the presence of a non-turbulent surface layer: vertical profiles of the radial (left) and azimuthal 
(right) velocity (upper panels) and magnetic field (lower panels) for the source term brs- Note that the radial velocity is normalised by 
aOi the diffusion coefficient inside the disc. The calculation using the physically motivated position of the non-turbulent layer (^cut = 4 
is shown with red lines, (^cut = 3 is shown with a blue line and the uniform diffusion coefficient case is shown in black. The vertical 
dotted lines show the height above which the magnetic pressure dominates over the thermal pressure. The magnetic field strength 
is /3o = 10*. The non-turbulent surface layer at its expected position barely cha nges the vertical profiles and the diff usion velocity. A 
significant bending of the field line inside the layer (as assumed in the model of iBisnovatvi-Kogan fc Lovelac3 l l2007t) 'l and associated 
reduction of the diffusion occurs only if the non-turbulent layer is artificially extended below its expected height (^b- 



tained only when the non-turbulent layer is artificially ex- 
tended to a smaller height, as shown with a blue line for 
Ccut = 3. In this case a significant part of the bending of 
the field line occurs in the non-turbulent layer (upper left 
panel). As a result the difi'usion velocity of the magnetic 
flux (visible as the radial velocity at large Q normalised 
by ao is decreased compared to the case of uniform dif- 
fusion coefficients. This result is further demonstrated by 
Figure 1131 showing the diffusion velocity of the magnetic 
flux normalised by ao as a function of the height of the 
non-turbulent layer (cut normalised by its expected value 
(b- The diffusion velocity of the magnetic flux is unchanged 
as long as C^cut > Cs. It is decreased by the presence of the 
non-turbulent layer only if C^cut < Cs • This is consistent with 
the averaging procedure described in Section [2]3l we argued 
that the average resistivity should be computed up to the 
height ("b, therefore the non-turbulent surface layer should 
extend below ("b in order to decrease the average resistivity 
determining the diffusion of the magnetic fiux. 

The numerical results presented above disagree with the 
expectations of [Si snovatv i-Kogan & Lovelace (20o3) and 
iRothstein fc Lovelace! (|2008l ). Indeed, we find a stringent 
constraint on the position of the non-turbulent layer in or- 
der to decrease the magnetic flux diffusion. This constraint 
is intimately linked with the condition on the magnetic field 
strength that should be fulfilled to enable a bending of the 
mag netic field lines. In Bisnovatvi-Kogan fc Lovelacd (|2007l ) 



and IRothstein fc Lovelacd 1 20081 ') this condition was evalu- 
ated by comparing the radial magnetic force to the grav- 
ity. This may however be a very optimistic choice because 



it would correspond to a non-rotating surface layer, which 
seems a bit extreme. The asymptotic expansion shows that 
the deviation from Keplerian rotation is of the order of 
Sv,jy ~ ■f^Cs ~ ^^vk- As a result of this smaller devia- 
tion from Keplerian rotation, the condition for a significant 
bending of the magnetic field lines on a length-scale of order 
_ff is: /3 > 1 or equivalently C ^ Cs- This is consistent with 
the constraint on the position of the non-turbulent layer 
Ccut < Cs so that the field lines can bend there. 

Note that with Ccut = Csi the above argument would 
still leave some hope that a significant bending could oc- 
cur in the non-turbulent layer if this bending can occur for 
/3 ~ 1. And the jump in br observed to occur around Cs 
shows that such bending can indeed occur. However, this 
bending is no longer controlled only by the diffusion as in 
the passive field region, but also by dynamical constraints 
since the magnetic force is no longer negligible. As a result it 
is crucial to self-consistently compute the velocity and mag- 
netic field profile to determine what bending can occur in 
this region. The numerical results show that the bending of 
the field line around Cs is always in the wrong direction: 
the radial magnetic field decreases around Cs after reach- 
ing a maximum at C ^ Cs- This feature is robust and was 
obtained in Paper I as well as Section |3] and |4l It was in- 
terpreted in Paper I as due to the fact that the azimuthal 
velocity at C ^ Cs is super-Keplerian when the magnetic 
field lines are bending outward {bra > 0). To enable super- 
Keplerian velocities the magnetic force should be directed 
inwards and therefore the radial magnetic field should de- 
crease with height. As a result (as shown in Figure [T2t when 
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Ccut — Cb, the non-turbulent surface is inefficient at reducing 
the diffusion. 



6 CONCLUSIONS 
6.1 Summary 

We have studied the dependence of the transport of mag- 
netic flux on the vertical profile of the diffusion coefflcients. 
For this purpose we have examined three different physically 
motivated vertical profiles. The first profile (Section[3|) is in- 
spired by the results of numerical simulations of MHD tur- 
bulent stratified discs. It therefore aims at describing fully 
turbulent discs in a more realistic way than the uniform dif- 
fusion coefficient model presented in Paper I. We confirm 
the main results of Paper I with this more realistic treat- 
ment, thus showing their robustness: the magnetic flux is 
advected more easily for weak magnetic field strength, re- 
sulting in a steady-state radial inclination of the magnetic 
field lines that c an be up to 30 times larger than predicted by 
previous works (|Lubow et al.|[l993 : iHevvaerts et al.lll99d ). 

In Section |4l we studied a model of disc containing a 
non-turbulent dead zone. We found that the presence of the 
dead zone does not change the radial inclination of the field 
lines in a steady state. It does change, however, the ratio 
of the advection velocity of magnetic flux to that of mass: 
the magnetic flux can be advected up to 100 times faster 
than mass for the dead-zone model described in this pa- 
per, and even more for more extended dead zones. This be- 
haviour could lead to an interesting time-dependence of the 
magnetic field configuration, which has the ability to evolve 
much faster than the mass distribution. 

In Section O we studied the effect of a non-turbulent 
surface layer on the diffusion of the magnetic flux. We found 
that it is able to reduce the diffusion of the magnetic flux 
only if it extends below the height (^b at which the magnetic 
pressure of the large scale vertical field equals the thermal 
pressure. Since the magnetic field is able to stabilise the MRI 
only above (b, the non-turbulent surface layer is inefficient 
at reducing the diffusion. 

All these results can be explained using the vertical av- 
eraging procedure described in Section [2]3l The relevant av- 
erage resistivity determining the diffusion of the magnetic 
field is the inverse of the volume average conductivity per- 
formed up to the height (b (where the magnetic pressure 
equals the thermal pressure). The averaged radial velocity 
relevant to the advection of magnetic flux is conductivity 
weighted (also performed up to the height ("s) and can there- 
fore be very different from the dens ity-weighted average t hat 
determines the transport of mass (|Ogilvie fc Lividl200ll ). 



6.2 Status of the 'magnetic flux problem' and 
future directions 

Early theoretical calculations predicted that the magnetic 
field diffuses faster than it is advected by the accre- 
tion flow, implying that a significant m agnetic flux should 
not be present i n thin accretion discs (jLubow et al.l[l993 : 
iHevvaerts eFall Il996l l. On the other hand, observations 
show the existence of strong outflows in a variety of objects 



containing an accretion disc (AGN, XRB, cataclysmic vari- 
ables, T-Tauri stars). Theoretical models explaining these 
outflows generally need a strong net vertical magnetic field 
to accelerate these outflows through th e magnetorotational 
mechanism l|Blandford fc Pavnd Il982l ) or the Blandford- 
Znajek mechanism for sy stems containing a black hole 
l|Blandford fc Znaiekll 19771 ). How could such strong vertical 
magnetic flelds be obtained? 

It has been proposed that the presence of a non- 
turbulent layer at the surface of the disc could solve 
the problem by drastically r educing the diffusion of the 
magnetic fleld jBisnovatvi-Ko gan fc Lovelace! |2007| . |2012| : 
iRothstein fc Lovelacell2008l : iLovelace et aLlbOOg T We have 
studied this possibility in Section [5] and found that the non- 
turbulent surface layer is inefficient at reducing the diffusion 
if it lies where the large scale magnetic fleld is dominant (as 
it is expected to). This conclusion should ideally be verified 
using numerical s imulations of turbulent discs. 

In Paper I l|Guilet fc Ogilvid lioi^ l. we showed that 
the situation was not that bad for weak magnetic fields 
iPo 3> 1). Indeed, the vertical structure of the velocity and 
magnetic field induce a faster advection and a slower diffu- 
sion of the magnetic flux as compared to the vertical aver- 
aging used before bv lLubow et al.l (| 19941 ) and others. As a 
consequence the inclination of the magnetic field lines in a 
stationary situation, which is regulated by the relative effi- 
ciency of advection and diffusion processes, could be 10 to 50 
times larger than thought before, for weak enough magnetic 
flelds. This encouraging result was obtained in the idealised 
case of uniform diffusion coefficients. In this paper, we have 
confirmed the robustness of this result with different vertical 
profiles of the diffusion coefficients, modelling either a fully 
turbulent disc or one containing a non-turbulent dead zone 
near its midplane. 

Though encouraging, these new results are so far in- 
sufficient for very thin discs, and obtaining strong magnetic 
fields here remains problematic. Additional effects neglected 
in the present study could help further the advection of the 
magnetic flux: 

• A more extended corona than the present isothermal 
model (either because it is hot or because of magnetic sup- 
port) would further reduce the diffusion (by increasing Cs for 
a given magnetic fleld strength) and possibly also increase 
the advection efficiency. 

• An anisotropic diffusion could also contribute to reduc- 
ing the diffusion of magnetic flux: if the vertical diffusion 
of the magnetic field is slower than the radial diffusion of 
angular momentum by some numerical factor, this would 
favour the advection of magnetic flux over diffusion. Such an 
anisotropy of the diffusion processes is suggested by the nu- 
merical results of iLesur fc Longarettil l|2009l ). though more 
numerical studies are certainly necessary to confirm this (the 
vertical diffusion of a radial magnetic field has never been 
measured so far, since iLesur fc Longarettil l|2009l ) only mea- 
sured the vertical diffusion of an azimuthal magnetic field). 

• Also we note that the diffusion properties of the mag- 
netically dominated corona have never been investigated. 
The importance of this region of the disc for the transport 
of magnetic flux ha s been highlighted by t his paper and Pa- 
per I as well as bv lBeckwith et al] l|2009l ). If this region of 
the disc was not efficient at diffusing magnetic field (though 
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it is highly dynamic), this could help further reduce the dif- 
fusion. 



Even with these additional effects something more may 
be needed to achieve strong magnetic field that can launch 
a strong outflow. In Paper I, we proposed that this ad- 
ditional process could be the outflow itself if some out- 
flow can be obtained for rather weak magnetic fleld fields. 
In this scenario the advection due to the effective turbu- 
lent viscosity would be efficient up to a certain strength 
of order (3 ~ lO'^ — 10* (depending on the aspect ratio of 
the disc). A (weak) outflow would then be launched and 
would enable further advection of the magnetic flux, flnally 
reaching strong magnetic flelds. The effect of an outflow 
on the ad vection of magne tic flux has been co nsidered by 
Lovelace _ et al. (2 009)j^I3is novatvi-Kogan fc Love lace (20l3)) 
Gullet fc Orilviel (|2012l ). A few numerical simulations sug- 
gest the possibility of such ou tflows with a weak mag- 
netic field. iMurphv et all (|201(]| ') have obtained a magneto- 
centrifugal outflow from a weakly magnetised disc with /3 ~ 
100 — 1000. An outflow was also obtained in the local shear- 
ing box model for even weaker m agneti c fields jP = 10"* — 
10^) b y ISuzuki fc Inutsukal ()2009( ) and lOkuzumi fc Hirosd 
l|201ll ). Much more work is, however, needed to understand 
the properties of outfiows that may be launched from weakly 
magnetised discs and the effect of such outflows on the trans- 
port of magnetic flux. 

We found in Paper I that the aforementioned scenario 
to obtain strong magnetic fields by inward flux transport 
could operate in discs with an aspect ratio above a (rather 
uncertain) critical value. This criterion may be met in proto- 
planetary discs, the aspect ratio of which is not very small: 
H/r ~ 0.05 — 0.1. Strong large-scale poloidal magnetic fields 
may then be obtained, enabling the launching of a magneto- 
centrifugal outfiow which could explain the jets observed in 
T-Tauri stars as argued bv lFerreira et al] ()2006h . In the soft 
state of X-ray binaries, the accretion disc is thinner (unless 
the luminosity approaches the Eddington luminosity), and 
is therefore likely to be too thin for the above scenario to op- 
erate. The disc would then be threaded by a weak magnetic 
field, which is not able to drive a powerful jet. This may 
explain the absence of jets in the soft st ate of X-ray binaries 
l|Fender et al.lll999l : ICoriat et aLll2011^ . In the hard state 
on the other hand, the inner part of the accretion flow is 
thought to be radiatively inefficient and geometrically thick. 
In principle, a strong magnetic field could therefore accumu- 
late in the inner geometrically t hick disc, and drive a power- 
ful jet such as simulated e.g. bv lTchekhovskov et al.l j201ll l. 
In order to reach this inner region, however, the magnetic 
fiux also needs to be transported inwards in the outer parts 
of the disc, which may be geometrically thin. The results 
of this paper and Paper I show that this fiux transport is 
possible for weak enough magnetic fields even if the disc 
is thin. One could therefore envision a situation where the 
outer thin accretion is weakly magnetised but allows mag- 
netic fiux to be transported inwards. This magnetic fiux can 
then accumulate and become dynamically significant in the 
inner geometrically thick accretion fiow, where a strong jet is 
launched. This idea should be tested in a global disc model, 
in which the time-evolution of the magnetic fiux could be 
computed using the transport rates computed in this paper. 



ACKNOWLEDGMENTS 

This research was supported by STFC. We thank Henk 
Spruit, Henrik Latter, Clement Baruteau, and Sebastien 
Fromang for useful discussions. 



REFERENCES 

Balbus S. A., Hawley J. F., 1991, ApJ, 376, 214 
Beckwith K., Hawley J. F., Krolik J. H., 2009, ApJ, 707, 
428 

Bisnovatyi-Kogan G. S., Lovelace R. V. E., 2007, ApJL, 
667, L167 

Bisnovatyi-Kogan G. S., Lovelace R. V. E., 2012, ApJ, 750, 
109 

Blandford R. D., Payne D. G., 1982, MNRAS, 199, 883 
Blandford R. D., Znajek R. L., 1977, MNRAS, 179, 433 
Coriat M., Corbel S., Prat L., Miller- Jones J. C. A., Cseh 
D., Tzioumis A. K., Brocksopp C, Rodriguez J., Fender 
R. P., Sivakoff G. R., 2011, MNRAS, 414, 677 
Fender R., Corbel S., Tzioumis T., Mclntyre V., Campbell- 
Wilson D., Nowak M., Sood R., Hunstead R., Harmon A., 
Durouchoux P., Heindl W., 1999, ApJL, 519, L165 
Ferreira J., Dougados C, Cabrit S., 2006, AfcA, 453, 785 
Fleming T., Stone J. M., 2003, ApJ, 585, 908 
Fromang S., Hennebelle P., Teyssier R., 2006, AfcA, 457, 
371 

Fromang S., Nelson R. P., 2006, A&A, 457, 343 
Fromang S., Terquem C, Balbus S. A., 2002, MNRAS, 329, 
18 

Fromang S., Terquem C, Nelson R. P., 2005, MNRAS, 363, 
943 

Gammie C. F., 1996, ApJ, 457, 355 

Gullet J., Ogilvie G. I., 2012, MNRAS, 424, 2097 

Heyvaerts J., Priest E. R., Bardou A., 1996, ApJ, 473, 403 

Kley W., Lin D. N. C, 1992, ApJ, 397, 600 

Lesur G., Longaretti P.-Y., 2009, AfcA, 504, 309 

Lovelace R. V. E., Rothstein D. M., Bisnovatyi-Kogan 

G. S., 2009, ApJ, 701, 885 
Lubow S. H., Papaloizou J. C. B., Pringle J. E., 1994, MN- 
RAS, 267, 235 
Miller K. A., Stone J. M., 2000, ApJ, 534, 398 
Murphy G. C, Ferreira J., Zanni C, 2010, AfcA, 512, A82 
Ogilvie G. I., Livio M., 2001, ApJ, 553, 158 
Okuzumi S., Hirose S., 2011, ApJ, 742, 65 
Rothstein D. M., Lovelace R. V. E., 2008, ApJ, 677, 1221 
Spruit H. C, 2010, in T. Belloni ed.. Lecture Notes in 
Physics, Berlin Springer Verlag Vol. 794 of Lecture Notes 
in Physics, Berlin Springer Verlag, Theory of Magnetically 
Powered Jets. p. 233 
Suzuki T. K., Inutsuka S.-i., 2009, ApJL, 691, L49 
Takeuchi T., Lin D. N. C, 2002, ApJ, 581, 1344 
Tchekhovskoy A., Narayan R., McKinney J. C, 2011, MN- 
RAS, 418, L79 
Terquem C. E. J. M. L. J., 2003, MNRAS, 341, 1157 

This paper has been typeset from a T^jiX/ I^T^tX file prepared 
by the author. 



